BEILSTEIN JOURNAL OF NANOTECHNOLOGY 



Sublattice asymmetry of impurity doping in graphene: 
A review 

James A. Lawlor*'' and Mauro S. Ferreira'''^ 



Review 






Address: 


Beilstein J. Nanotechnol. 2014, 5, 1210-1217. 


''School of Physics, Trinity College Dublin, Dublin 2, Ireland and 


doi:10.3762/bjnano.5.133 


^CRANN, Trinity College Dublin, Dublin 2, Ireland 






Received: 29 April 2014 


Email: 


Accepted: 08 July 2014 


James A. Lawlor* - jalawlor@tcd.ie 


Published: 05 August 2014 


* Corresponding author 


This article is part of the Thematic Series "Atomic scale interface design 




and characterisation: Theory - Electronic properties". 


Keywords: 




band gap; electronic transport; graphene; nitrogen doping; sublattice 


Guest Editor: A. Krasheninnikov 


asymmetry 






©2014 Lawlor and Ferreira; licensee Beilstein-lnstitut. 




License and terms: see end of document. 



Abstract 

In this review we highlight recent theoretical and experimental work on sublattice asymmetric doping of impurities in graphene, 
with a focus on substitutional nitrogen dopants. It is well known that one current limitation of graphene in regards to its use in elec- 
tronics is that in its ordinary state it exhibits no band gap. By doping one of its two sublattices preferentially it is possible to not 
only open such a gap, which can furthermore be tuned through control of the dopant concentration, but in theory produce quasi- 
ballistic transport of electrons in the undoped sublattice, both important qualities for any graphene device to be used competetively 
in future technology. We outline current experimental techniques for synthesis of such graphene monolayers and detail theoretical 
efforts to explain the mechanisms responsible for the effect, before suggesting future research directions in this nascent field. 



Review 
Introduction 

With its excellent transport properties and low dimensionality, 
graphene, an atomically thin layer of Carbon atoms bonded 
together in a hexagonal lattice, initially seems a strong candi- 
date for use in many future commercial applications such as 
ultra high-speed transistors, integrated circuits and other novel 
devices [1-3]. One of the main problems with using regular 
graphene for such applications is the absence of a band gap in 
the electronic band structure [4], and as a result any field effect 



transistors (FETs) made using the material (so-called GFETs) 
would be unable to be switched off, rendering it useless as a 
logic device [5-7]. The most natural way to approach this issue 
is therefore to introduce a sizeable band gap and hence allowing 
more control over the current flow. 

Although several methods exist to induce a band gap, for 
example ID quantum confinement by construction of graphene 
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nanoribbons (GNRs) [8], stacking of monolayers with perpen- 
dicular electric fields [6,7], strain [9] and mounting on a sub- 
strate [10-12], these methods are not without problems. Alter- 
natives are therefore sought after with the minimum standard to 
meet or exceed the limits of silicon semiconductor technology, 
which is characterised by a current on/off ratio of roughly 
10^-10^ and a band gap of at least 340 meV whilst maintaining 
high carrier mobility [5,13,14]. 

Alteration of the crystal structure through the introduction of 
foreign dopants is one of the more realistic avenues of ap- 
proach in realising this goal. Atomic dopants like boron (B) or 
nitrogen (N) are a similar size to carbon and can be introduced 
easily in a variety of graphene growth processes, typically 
replacing carbon sites and forming substitutional impurities in 
the lattice, positively (p-) and negatively (n-) doping the system 
for B and N dopants respectively [15]. Early theoretical 
attempts at investigating the electronic properties of such a ma- 
terial found that a periodic arrangement of B or N dopants, 
forming a dopant superlattice, would open a band gap [16], but 
that a random distribution of dopants among lattice sites yields 
no band gap [17]. Further investigations using vacancies, where 
carbon atoms are removed from the lattice, found that both 
superlattices [ 1 8] and random distributions restricted to one of 
the two graphene sublattices [19] both lead to a tunable band 
gap, and in the latter case an emergence of magnetic properties 
in the system [20,21]. 

Whilst these findings are certainly interesting, their scalability, 
and hence commercial application, is prohibited by the stan- 
dard of precision that must be met in order to produce such ma- 
terials. Circumventing these problems using nitrogen doped 
graphene has only come about relatively recently which we 
detail in the following section. 

A brief history of nitrogen doped graphene 

Exploiting the effects of nitrogen dopants on the transport prop- 
erties of graphene has been an interesting experimental research 
topic for the last 5 years [22], whereby the methods and tech- 
niques to achieve this have been developed to include chemical 
vapour deposition (CVD), using NH3 as a precursor, arc 
discharge [23], embedded nitrogen and carbon sources within a 
metal substrate [24], ion implantation [25,26], ammonia [27] or 
nitrogen plasma [28,29] treatments, and capable of achieving 
nitrogen dopant concentrations of up to around 10% [30] and 
with direct applicablity to GFET technology and bio-sensing 
[28]. Although CVD is one of the more challenging methods, it 
seems the most reliable option and yields the best quality 
nitrogen doped graphene sheets [31] and single continuous 
sheets can be synthesised on the centimeter scale [32]. Using 
nitrogen dopants alone through CVD can yield bandgaps up to 



200 meV [33] and by inclusion of boron co-dopants, through a 
tailored growth process, this can be expanded to around 
600 meV with 6% total dopant concentration [12]. Boron 
doping alone has shown to also open a band gap, tunable with 
dopant concentration [34,35]. This method results in small B-N 
domains embedded in the graphene sheet, so the precise mecha- 
nism behind the band gap opening may not be the same as that 
seen in a random B-N ensemble. Although these gaps are quite 
sizeable, the effect of scattering would lead to a detrimental 
impact on the transport characteristics. However, the recent 
theoretical prediction of sublattice asymmetry, a situation where 
dopants are distributed in one sublattice only, is expected to 
yield very low scattering transport properties. The electron 
wavefunctions in such a system are shown to mainly exist on 
the sublattice without dopants and can travel almost unhindered 
[36,37], which is very promising for overcoming the scattering 
problem. 

It was not until 201 1 that a viable experimental approach was 
found by Zhao et al. [38]. They discovered that graphene grown 
via chemical vapour deposition (CVD) in the presence of 
ammonia (NH3) naturally incorporates nitrogen atoms as substi- 
tutional so-called 'graphitic' dopants (see Figure lA) into the 
crystal, and with a distinct sublattice segregation of dopants. 
Indeed further research has uncovered a less pronounced asym- 
metry phenomenon using graphene implanted with nitrogen 
impurities followed by a high temperature annealing process 
[25], and it seems reasonable that there is a common mecha- 
nism with the CVD method. Due to the limited published work 
on this implantation and annealing method, our main focus will 
be on CVD and it is here that we begin our main discussions. 
We begin with a section covering in-depth the experimental 
work done to date, followed by a section on theoretical aspects 
and ending with a future outlook and conclusions. 




Figure 1 : Schematic of a graphene lattice with the most common 
experimentally observed species of substitutional nitrogen dopants: 
(A) a single graphitic, (B) three pyridinics, (C) one N2'^ pair, (D) one 
N2''^^ pair and (E) one N2^^ pair. 



1211 



Beilstein J. Nanotechnol. 2014, 5, 1210-1217. 



Experimental observation of sublattice asym- 
metry of N-graphene 

The main focus of the work by Zhao et al. was the observation 

and characterization of nitrogen dopants via scanning tunnel 
spectroscopy (STS), followed by a short investigation on the 
transport properties of the resulting graphene. Their experi- 
mental observation of same sublattice segregation, at least on 
local scales, was noted as a curiosity but was not discussed in 
depth. A key finding by the researchers was that by varying the 
ammonia precursor concentration they could adjust the resulting 
embedded dopant concentration. Further theoretical work 
confirmed that a tunable band gap would be obtainable using 
this method [36], and along with the findings of Zhao et al. this 
sparked further interest into this phenomenon. It should be 
noted that the experimental conditions were a high vacuum 
CVD process on a Cu(lll) substrate with a temperature of 
1000 °C, using CH4 and NH3 as graphene and nitrogen precur- 
sors respectively for 1 8 min reaction time achieving a dopant 
concentration of 0.3%. 

Shortly after this work a more focussed investigation was 
undertaken by Lv et al. [39]. They were able to replicate the 
findings of the first group and synthesised large areas of 
N-doped monolayer graphene with the additional finding that 
by increasing the pressure in the growth phase from a high 
vacuum to atmospheric-pressure CVD would result in an abun- 
dance of so-called N2'^'^ pairs [39] (see Figure IC), which are 
thought to be a result of a higher number of intermolecule colli- 
sions during formation. Using the same substrate and precur- 
sors, they were also able to find that a minimum reaction time 
of 5 min and temperatures of 800 °C were optimum for the 
N-doped graphene synthesis and moreover it was found that 
singular nitrogen dopants are more frequently found for reac- 
tion times below 10 min with the ratio of N2'^'^ to single N 
increasing with reaction time. The samples had a typical 
concentration of around 0.25% nitrogen dopants. 

More recently the addition of boron dopants into the lattice and 
their effects have been studied [40] using high vacuum CVD 

growth and adding B2H6 diborane instead of ammonia as the 
dopant precursor, achieving concentrations on the order of 
0.3%. The researchers compared B-doped and N-doped systems 
in detail and found no detectable sublattice asymmetry in the 
case of B-doped systems and thus the dopants were distributed 
evenly between sublattices. The reason for this is thought to be 
due to a strong interaction between the Cu(lll) substrate and 
the boron impurity, in contrast to a very weak interaction with 
the substrate for nitrogen impurities [40]. 

The most recent and thorough study of asymmetric nitrogen 
dopants in graphene is by Zabet-Khosousi et al. [41] who again 



used the CVD growth process but with a pyridine (C5H5N) 
precursor for the nitrogen and graphene instead of the conven- 
tional ammonia/methane mix. Figure 2 shows a rather striking 
STM image of a large area of N-doped graphene grown using 
this method, showing two clearly defined domains of different 
sublattice preference where the sharp borders are thought to 
arise from terracing of the substrate disrupting the asymmetry 
effect. The method was tailored to use a very smooth Cu(l 11) 
substrate and so they were able to achieve domain sizes far 
larger than previous efforts. Typical concentrations found with 
this method were around 0.2%, very similar to the previous 
ammonia/methane method which points to an equivalent growth 
mechanism. By comparing N-doped graphene systems synthe- 
sised with this method to other systems made via nitrogen ion 
bombardement and ammonia post-treatment of a pristine 
graphene sheet, where no sublattice asymmetric configurations 
were found, it was suggested that the sublattice asymmetry must 
be occuring during formation of the graphene sheet. This is 
supported by the recent experiments using ion bombardement 
followed by high temperature annealing by Telychko et al. [25], 
where the impurities may be reconfiguring themselves locally 
during cooling. This phenomenon is not seen without the 
annealing process [26]. 

The common elements of all the aforementioned experimental 
reports are that they all use a Cu(l 1 1) substrate, all have reac- 
tion times exceeding 10 min in a temperature range of 800-900 
°C, and all have used either ammonia or pyridine precursors to 
achieve concentrations of up to 0.3%. A review of other 
methods for N-doped graphene synthesis finds no mention of 
sublattice asymmetry [29,30], it is known that this effect is 
mainly detectable through careful scanning tunnel microscopy 
of large areas. The most natural explanation for this lack of 
observation is then either that it is not being observed due to 
simply not looking, or that sublattice asymmetry is not a very 
robust effect and therefore requires a careful choice of syn- 
thesis procedure. Research into the mechanisms behind the 
phenomenon is therefore paramount in answering this question 
and could shed light on whether other methods or even dopant 
species are at all possible and this leads us naturally to a review 
of the current theoretical models of the effect. 

Theoretical models of the segregation effect 

It is clear that any degree of asymmetry between the two 
equivalent sublattices is the result of a symmetry breaking oper- 
ation. Current theoretical attempts to explain the sublattice 
asymmetry in the nitrogen doping seen in the experiments 
discussed in the previous section suggest that the effect comes 
from either the energetically preferable positioning of nitrogen 
on the graphene edge during the growth process [41,42], where 
the symmetry breaking effect is the edge structure, or from 
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Figure 2: STM images of nitrogen doped graphene on (a) 7 nm^, (b) 20 nm^ and (c) 100 nm^ scales, adapted with permission from Zabet-Khosousi 
et al. [41]. Copyright 2014 American Chemical Society. The dopants locations are identified by finding bright spots on the STM image, corresponding 
to slight perturbations in the positions of the neighbouring carbon atoms. The dopant sublattice can be found through the orientation of the bright 
triangle features, where opposite sublattices appear as mirror images of each other as demonstrated in (a). The red and blue triangles in each 
subfigure correspond to impurities on the different sublattices. Part (c) best demonstrates that although the distribution of dopants appears random, 
when their sublattice is identified we see two very large domains appear with opposite sublattice segregation. Shown at the bottom, off-center left of 
(c) is a white 10 nm scalebar 



inter-impurity interactions in the impurity ensemble [39,43] 
where the symmetry is broken by the impurities themselves. 
These models will hereafter be referred to as the edge growth 
and interaction models respectively. 

The edge growth model suggests that during the growth phase 
the most energetically favourable position for a nitrogen dopant 
being incorporated into a graphene edge would be that which 
leads to a same sublattice configuration for all impurities in a 
domain. Through density functional theory (DFT) calculations 
involving a graphene nanoribbon on a Cu(lll) substrate, 
aiming to reproduce experimental conditions, a thorough 
investigation into the energetic favourable position of single 
graphitic nitrogens in the GNR was undertaken and the ener- 
getic difference between placing the nitrogen at one edge site 
over another on the opposite sublattice was found to be a 
substantial 1.3 eV. One drawback of this method is that one 
would expect therefore that the clearly defined segregation 
domains seen experimentally (see Figure 2 and [41]) would 
naturally fall in with the graphene grain boundaries, but this has 



not been observed. The role of inter-impurity interactions has 
been considered as an alternative, through both tight binding 
[43] and DFT [39,44] formalisms, although with differing 
conclusions. The tight binding method by Lawlor et al. [43] was 
part of a more in-depth theoretical investigation, and suggests 
that a Friedel oscillation-like perturbation in the long-range 
inter- impurity interactions [45] arise from adding impurities to 
the system, leading to a shift in the system's Fermi energy. 
Such a mechanism would explain both the non-commensura- 
bility between crystal grain boundaries and segregation 
domains, and would suggest why the alternative synthesis 
method of nitrogen implantation followed by high temperature 
annealing also results in impurities preferring to occupy the 
same sublattice [25]. Furthermore, the authors predict that as a 
result of this interaction and the reduction in nearest neighbour 
distance with increasing concentration, there exists a critical 
dopant concentration beyond which no sublattice asymmetry 
would be observed, as the energetic minimum occurs when the 
impurities are distributed evenly between sublattices. Although 
the exact value of this critical concentration is dependant on 
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how the impurity is parameterised within the tight binding 
regime, the authors predict it to lie between 0.1% and 0.8%. 
Comparing to the highest experimentally reported doping 
concentrations for samples with sublattice asymmetry (0.3%), it 
is a hope that further experiments can clarify whether such 
a critical concentration exists. If the prediction of a critical 
concentration is accurate it would limit the band gap of a 
segregated device to around 100 meV [36]. This figure, 
however, comes from matching tight binding and DFT band 
structure results, a method which is known to systematically 
underestimate such band gaps. Nevertheless, the band gap 
obtained can be expected to be much below that required for a 
GFET device. 

In-depth DFT calculations by Hou et al. [44] found that the 
interactions between nitrogen impurities are generally repulsive 
in nature. More specifically, when the two impurities are placed 
close together the system energy is minimised for opposite 
sublattice configurations, with the exception of N2'^^. This 
contradicts the experimental reports of Lv et al. [39] who found 
an abundance of N2^^ defects in their samples. A simplistic 
tight binding approach predicts that two identical impurities in 
close proximity to one another would indeed have an energetic 
minimum when both occupy the same sublattice [46], however 
this method ignores coulombic interaction which will be domi- 
nant at close separations. Interestingly, the DFT approach of Lv 
et al. of such nitrogen pairs finds that the lowest energy con- 
figuration would be N2'^^ (see Figure 1), approximately 0.3 eV 



lower than that of N2 . This appears to be in disagreement 
with experiments where N2'^^ is more commonly found than 
N2^^ . As a final note their calculations have also shown that 
two N2^^ pairs have a lower energy when they share the same 
sublattice, however further studies have shown that the overall 
energy change resulting from an N2^ pair in the sheet is very 
high [42]. It is possible that these impurities also come from 
edge growth, but the controllable presence of N2^'^ in nitrogen 
ion implanted graphene [25] indicates again that this may not be 
the full picture as in this case the complete graphene sheet is 
already fabricated. 

Predicted electronic properties 

The earliest attempts to study nitrogen dopants and their effect 
on the graphene electronic properties were purely theoretical. 
Beginning with impurity superlattices, this research preceded 
the experimental realistation of nitrogen doped graphene. By 
introducing a controlled periodic arrangement of boron and 
nitrogen impurities a band gap was seen to open [16], and 
although such superlattice structures were not feasible on large- 
scale it was further shown that a random distribution on one 
sublattice can also open a gap [37]. Further DFT studies showed 
that the band gap will increase with dopant concentration 
[36,47] and that a dopant level of over 8%, where impurities are 
all on the same sublattice, will produce a band gap of around 
550 meV far surpassing the minimum required for a CMOS 
[14,36] and finding that the band gap scales with concentration 
to the power 3/4, as shown in Figure 3. 




Figure 3: Predicted band gap against nitrogen dopant concentration. Blacl< circles are calculated values from [36] and tine red dashed line shows the 
expected band gap scaling with concentration, according to the power 3/4 as discussed in the text. 
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Even with a 4: 1 doping ratio between sublattices, the band gap, 
although smaller, was shown to still exist. This is promising for 
scaleability where perfect asymmetric doping may not always 
be realizable. Beyond the realization of quasi-ballistic electron 
transport, sublattice segregated systems can also be used to in- 
duce magnetism [19,48-50] and produce spin-polarized current 
[51,52]. 

Much work has been done studying how the placement of 
dopants affects the properties of nanoribbons. DFT approaches, 
using a periodic system of dopants [53,54], and a more general 
Kubo-Greenwood approach [37] have shown that electron 
transport is enhanced when dopants are placed on one sublat- 
tice, compared to a random distribution, and that a band gap 
does indeed open. The difference in transport qualities between 
the asymmetrically doped versus completely randomly doped 
systems is illustrated by a conductance plot in Figure 4, where it 
is evident that the electrons in the fornier kind of device will be 
subjected to less scattering leading to an increase in the 
quantum conductance, closer to that of the pristine system when 
subjected to a positive bias. 

It should also be noted that that symmetry breaking in nanorib- 
bons occurs via edge effects and also results in a favourable 
sublattice [58], however the mechanism is distinctly different 
from that in graphene. 



Outlook 

At the time of writing, the synthesis of sublattice asymmetric 
graphene is only just now becoming practical. There are still 
many unanswered questions and the most important ones we 
outline in this section. 

The natural question to ask is if this asymmetry is particular to 
nitrogen only, or if it can be found with other dopants. Theoreti- 
cal findings in the past have shown that such an effect can be 
expected with a dilute concentration of certain adsorbates [59], 
and indeed more recent studies suggest that this could also be 
possible with other impurities [43]. Whilst it is known experi- 
mentally that molybdenum impurities exhibit same sublattice 
configurations in bilayer epitaxial graphene, the mechanism 
behind this is not currently understood [60]. Boron has been 
studied in a similar CVD growing regime to nitrogen where it 
has been proposed that the role of the substrate can play a crit- 
ical role in the manifestation of segregation. In this case, it is 
thought that the strong interactions between the Cu(l 1 1) crystal 
and the boron dopants destroy any asymmetry effects [40]. It is 
then logical to ask whether suitable substrates can be identified 
to produce segregated boron doped graphene sheets, which 
would have the effect of p-doping the system, and then to ask 
whether this can be extended to other species of dopants. We 
note that it is known that there are weak interactions between 
graphene and Al, Ag, Au, Pt(l 1 1) substrates, all of which leave 




Figure 4: Quantum conductance through a 1 5 nm wide graphene nanoribbon with a 7.5 nm long scattering region containing a dispersion of substitu- 
tional nitrogen impurities, in a similar vein to the method of Botello-Mendez et al. [37], calculated using a recursive Green's Function method [55,56], 
the Kubo formula for conductance [57] and a configurational average of 50 systems. Energy is in units of the tight binding nearest neighbour hopping 
energy between carbon atoms, f = 2.7 eV. Shown is the predicted conductance for pristine (black), randomly doped (green solid) and single sublat- 
tice doped (red dashed) systems, where a dopant concentration of 1% nitrogen was used. 
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the electronic structure intact [61], whilst substrates such as Ni 
have considerably stronger interactions [62]. Another way to 
investigate the presence of sublattice asymmetry with dopants 
other than nitrogen is via ion bombardement [63], which could 
be combined with the high temperature annealing process 
discussed previously. 

The experimental realisation of spin-polarized transport should 
also be pursued, and could spark additional interest in this 
research area beyond quasi-ballistic transport. Investigation of 
the effects of strain on such a system is difficult due to the tech- 
nicalities of the CVD method, but could be explored using 
adsorbates instead of substitutions or by incorporation of the 
strain in to the nitrogen ion bombardement and annealing proce- 
dure discussed earlier. Another open question is if a critical 
concentration of dopants exists, as mentioned in the theory 
section, which would limit the available bandgap below the 
threshold needed for future use in GFET devices. Currently all 
experimental reports on sublattice asymmetry have very low 
concentrations of dopants around 0.3%, so it should be feasible 
to test the existence of a critical concentration of using avail- 
able methods [30]. Consequentially this would shed more light 
on the mechanism responsible for the segregation, and whether 
in fact the inter- impurity interactions and edge growth effects 
are complementary effects. 

Conclusion 

This paper reviews the current state of experimental and theo- 
retical research in sublattice asymmetric nitrogen doped 
graphene. Not only are such systems now able to be synthe- 
sised in the lab, we have seen that this area shows promise 
attaining graphene-based FET devices due to the opening of a 
bandgap whilst maintaining the excellent transport properties of 
graphene, something which is not realisable in graphene with no 
sublattice imbalance in dopant distribution. While there are still 
many open questions in the field these should be answerable 
within the scope of current techniques. 
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